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We find a mapping of the layered sine-Gordon model to an equivalent gas of topological excitations 
and determine the long-range interaction potentials of the topological defects. This enables us 
to make a detailed comparison to the so-called layered vortex gas, which can be obtained from 
the layered Ginzburg-Landau model. The layered sine-Gordon model has been proposed in the 
literature as a candidate field-theoretical model for Josephson-coupled high-T c superconductors, 
and the implications of our analysis for the applicability of the layered sine-Gordon model to high- 
Tc superconductors are discussed. We are led to the conjecture that the layered sine-Gordon and 
the layered vortex gas models belong to different universality classes. The determination of the 
critical temperature of the layered sine-Gordon model is based on a renormalization-group analysis. 
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I. INTRODUCTION 



Two essential prerequisites for an analysis of superconducting materials are anisotropic models, as initiated by 
Ginzburg [l[, and the inclusion of vortices, as envisaged by Abrikosov [2J. Typical high transition temperature super- 
conductors consist of copper-oxide superconducting planes separated by insulating layers. In the phenomenological 
description of high-Tc superconductivity, one may use an anisotropic, continuous Ginzburg-Landau theory [3, [3,0], 
but only for not too large anisotropy. Note that the anisotropic, continuous model can be mapped onto the isotropic 
Ginzburg-Landau model by an appropriate rescaling method [5j. However, in the case of extremely high anisotropy 
like in Bi2Sr2CaCu20s, the discreteness of the structure becomes relevant jf|, and it becomes necessary to use a 
layered Ginzburg-Landau model 0, 0] where the layers are coupled by Josephson or electromagnetic interactions. 
This provides a good basis for the discussion of the vortex dominated properties of high-T c superconductors. Some 
exact and some approximate mappings of the layered Ginzburg-Landau model (i.e. the Lawrence-Doniach model Q) 
onto various other statistical, field theoretical or spin models, like the layered vortex gas [H, [H, [S ES El, E3 and the 
anisotropic XY models [TH, EH, EH have already been discussed in the literature, and these models have also been 
proposed and used for the description of the vortex dynamics in high-T c superconductors. Connections to sine-Gordon 
type models [H El ES IS HS HH, S3, H3 have also been explored in the literature. 

The latter investigations are motivated by the well-known fact that the massless two-dimensional (2D) sine-Gordon 
scalar field theory belongs to the universality class of the 2D-XY spin model and consequently to that of the 2D 
Coulomb or vortex gas. The mappings between these models and also the phase structure have been discussed in 
the literature in great detail (see, e.g., Refs. @, ©, ES El E3, El, El El, HI El, HS El HS HI])- Sincc the layered 
Ginzburg-Landau model can be considered as the continuum limit of the anisotropic 3D-XY model (discrete in the 
z-direction but continuous in the xy- planes) , one might suggest that the field theoretical counterpart of the layered 
Ginzburg-Landau model should be a sine-Gordon type model. However, the 3D Ginzburg-Landau theory (in the 
London limit and in the absence of electromagnetic fields), which can be considered as the continuum limit of the 3D- 
XY planar rotator, and the 3D sine-Gordon model do not belong to the same universality class (see Refs. [30ll3lll32]|). 
a phase transition being absent in the 3D sine-Gordon case. Since layered models are always constructed from 3D 
models by a suitable discretization of the derivative term in one of the spatial dimensions (see, e.g., Ref. [HI), the 
equivalence of the layered vortex gas and layered sine-Gordon models remains questionable. One purpose of this 
paper is to clarify this point by finding an exact mapping of the layered sine-Gordon model to an equivalent gas 
of topological excitations, which in turn can be compared directly to the layered vortex gas. We also consider the 
phase structure and the critical behaviour of the iV-layer sine-Gordon model by a renormalization group method, and 
determine the relation of the critical parameter 6 2 of the layered sine-Gordon model to the critical temperature, as a 
function of the number of layers. 

This paper is organized as follows. In Sec. |H] we discuss the comparison of the layered Ginzburg-Landau to the 
layered sine-Gordon model, by a mapping of each model to an equivalent gas of topological excitations. In Sec. IIII1 
we discuss the the renormalization-group (RG) flow of the layered sine-Gordon model. Conclusions are reserved for 
Sec. EH 
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II. LAYERED GINZBURG-LANDAU VERSUS LAYERED SINE-GORDON MODEL 
A. Mapping of the layered Ginzburg— Landau model to the layered vortex gas 

The Ginzburg-Landau theory has been developed by applying a variational method to an assumed expansion of 
the free energy in powers of \ip\ 2 and |<9 At '0| 2 where ip is a complex order parameter (the inhomogeneous condensate 
of the superconducting electron pairs) and \ip\ 2 represents the local density of superconducting electron pairs (for a 
detailed discussion see, e.g., Ref. [HI). Its detailed form can be found in Eqs. (6-6) and (6-9) of Ref. [35|. Upon a 
discretization of one of the spatial directions (say, the z-coordinate), one obtains the layered Ginzburg-Landau (or 
Lawrence-Doniach Q) model with the free energy (in natural units: h = c = eg = 1) , 



^rl^r 1 + — 2^ — J + e 



\4>n+l ~ Ipnl 2 

2 m r 



(1) 



Here, m a b and m c represent the intralayer and interlayer effective masses, s is the interlayer distance, and N stands 
for the total number of layers. The sum over /x covers the spatial coordinates fi = x,y,z. The parameters a and (3 
are discussed in Eqs. (6-6) and (6-8) of Ref. [351 ] - In order to investigate the vortex dynamics in the framework of the 
Ginzburg-Landau theory, one has to consider the discretized model given in Eq. |T]) in the London approximation. 
One writes the complex, layer-dependent order parameter as ij) n {r) = V'On(^) exp[i</>„(r)] with real ipo in (r), where 
the 4> n € [0, 27r) are compact variables, and the moduli ipo t n are assumed to be constant and identical in every layer 
(i.e. ipo,n('f f ) = V'o) which is the London-type approximation. The London-type form of the layered Ginzburg-Landau 
model with Josephson coupling can be mapped (see Refs. [1, d, [j| [l(| El, EH) onto the layered vortex-gas. The 
globally neutral layered vortex-gas with N layers is characterized by the partition function (see Eq. (2.3) of Ref. (Tlj |) 

N - j2 N 



Z ^g = J2 p)2 E J ■■■ J2J E exp I -— ^-a Q ^y(r^,n a ^) I , (2) 

v ^ ni 1 ri2 v 1 cri,...,fT w =il V a i^0 / 

where a a = ±1 is the charge of the ath vortex, a stands for the lattice spacing, k-Q is the Boltzmann constant, T is the 
temperature, and the interaction potential V between two vortices depends on their relative distance r a p within the 
two-dimensional planes (r a p = \ f a — rp\) and on the distance n a p across the planes (n a p = \n a — n^Q, where n a is the 
layer in which the ath vortex is located. There are 1v vortices with fugacity z and these fulfill the neutrality condition 
y] T a a = . The positive (negative) vorticity is represented by positive (negative) charges. The restriction to 
globally neutral char ge c onfigurations is ensured by the condition ov+ 7 = — er 7 for 7 G {1, ... v\ in Eq. ([2|). Following 
Refs. @, H, H, 0, [Hljll2|' we neglect interactions between vortices separated by more than one layer, and this results 
in intra- and interlayer interaction potentials which have commonly accepted short- and long-range asymptotic forms 
given by (see, e.g., Eqs. (2.5) and (2.6) of Ref. p]]) 

V{r aPl n a!i = 0) = - In r J*^±, (3a) 

V a / a 

V(r aP ,n aP = l) = 6VA^. (3b) 



The coupling A ~ a 2 J±/J\\ is proportional to the ratio of the interlayer Josephson coupling J± to the intralayer 
coupling J11 , and b is a constant of order unity. The intralayer interaction between the vortices is logarithmic for 
short distances, as in the case of the usual 2D Coulomb or vortex gas, but linear for large distances. The interlayer 
interaction is always linear and similar to the long-range intralayer interaction but with an opposite sign. Within a 
layer, vortices of opposite charge attract, whereas the positive prefactor of the linear term in the interlayer interaction 
implies the formation of vortex stacks of like charges. 

B. Mapping of the layered sine— Gordon model to an equivalent gas of topological excitations 

The well-known sine-Gordon model in Euclidean space is defined via the action 

(4) 



Ssg[<p] = J d 2 r Q {d^ip) 2 ~ycos(bip)^j 
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where the minus sign of the periodic term is chosen so that the zero-field configuration remains a (local, not global, 
infinitely degenerate) minimum. As usual, ip here is a dimensionless scalar field, {d^tp) 2 = (<9 M c^) 2 , y is a funda- 

mental Fourier amplitude, and b is a dimensionless frequency. This model is well known to describe the Kosterlitz- 
Thouless-Berezinskii (KTB) phase transition 36] in two dimensions. If one adds a second layer which leads to the 
appearance of two fields tpi and tp2, one may devise the following natural ansatz for the interlayer interaction term 
i J(ipi — P2) 2 , where J is the Josephson-type coupling whose physical dimension is equal to the square of the inverse 
leng th. Indeed, the layered sine-Gordon model with this particular interlayer interaction term has been proposed in 
[Tft for the description of the vortex properties of Josephson coupled layered superconductors. The double-layer 

sine-Gordon model [la, E3, ii [S3, m 

is thus given by the Euclidean action 
SaLSG = J d 2 r ^(d M <p„) 2 + ijfai - tp 2 ) 2 - yJ2cos(b(p n ^J 

= J d 2 r ^(^) T (^) + ^ T g V-j/^co S (6/^)j , 



(5) 



where (p denotes the column vector ip = (^1,^2) characterizing the 0(2)-doublet, and the / are projectors / = 
{Si n , <>2n) whose components are given by Kronecker-deltas. The Josephson-type interlayer interaction corresponds to 
the following dimensionful mass matrix (see, e.g., Ref. pp||). 

- = (-J (6) 

The mass-eigenvalues are and 2J. In order to perform the mapping of the double-layer sine-Gordon model ([5]) onto 
a gas of topological excitations, we follow the scenario of Rcf. [3GJ, where the partition function of the sine-Gordon 
model is identically rewritten in the form of the partition function of a Coulomb gas. We should perhaps note that 
this mapping procedure is inspired by the treatment in Chap. 31 of [26| . One expands the exponential factor of the 
integrand with the periodic potential in a Taylor series, expresses cos(6/ T (/?) in terms of exponential functions and 
introduces integer valued variables, the charges a n = ±1, that fulfill the neutrality condition. After these operations, 
one obtains 

2 2LSG = Af I [%]exp(-5 2LSG M) = Aff^^t- £ f d 2 n . . . £ f d 2 r 2 „ 
J „=o y lv >- ni =l J n 2v =l J 

E / M CX P (- / d * r {\f (-± W + m 2 )^ + i bp T tp) ) , (7) 

CTi ,...,(T W — ±1 
CT„ + 7 = -cr T , 7g{l,...I/} 

where 1 stands for the two-dimensional unit-matrix which will be suppressed in the following. The charge density p(r), 
which depends on the configuration of the charges u\ , . . . , 02v and on their positions r\, . . . , ?2v, constitutes a vector 
in the internal space of the fields (tpi,ip2) characterizing the two layers and reads, p(r) = ^2 a=1 o a S(r — r a ) f . 
We have thus obtained a representation in which the 2v charges have been placed onto the two layers, with the ath 
charge on layer n a . Performing the Gaussian path integral in Eq. ([7]), we obtain 



£T;L,...,OV = ±l 

jv+ 7 =-«t t , 7G{l,...i/} 
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where p(p) = X)q=i °e» cxp(ip • r a ) f is the Fourier transform of the 0(2)-charge density. In momentum space, the 
propagator can easily be calculated by matrix inversion, and this gives 



--2LSG 



v=0 v ' ni=X 1 J 



E 



n 2v = l ' 



(2d 2u ^ 
-b 2 ^ ^ ^ a a (T 7 ( £>i„ Q 
CK — 1 7—1 



(71 ,...,(Tj,— ±1 

<x„+ T = -cr T , 7£{l,...l/} 



<$2n 



-<4('"a7) 

B{r ai ) 



B{r al ) 



Here, the interaction potentials are (r Q7 = \ f a — r 7 |) 

d 2 p gpp'-Cfc-^)] (p>2 + J) 



B(r ai ) = 



(2tt) 2 p 2 (p 2 + 2J) 
d 2 p elip-^-^)] J 
(2tt) 2 p 2 (p 2 + 2J) = ' 2n \2 



1 
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' 07 
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K 
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Acft 



' c*7 
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a 

A ff 



(9) 

(10a) 
(10b) 



where the momentum integrals can be performed using either dimensional regularization [26[ or ultraviolet (UV) cutoffs 
and the physically relevant, finite parts of the interaction potentials consist of massless and massive scalar propagators. 
In the expression for the intralayer (A) and intcrlayer (B) interaction potentials, a is the lattice spacing which serves 
as a short-distance (UV) cutoff, Kq denotes the modified Bessel function of the second kind, and A e g = 1/V2J is an 
effective screening length. The asymptotics of the interaction potentials read as follows (js = 0.55721 56649 ... is 
Euler's constant), 



A{r ai < Acff) 

A{r ai > A c ff) 
B{r ai < A cf t) 
B(r ai > A cf t) 



-Am 

2tt 



a 



0. 
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1 

27 
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A ff 
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ln(2) 



^ln(2)-i 7E 



7e 



(11a) 

(lib) 
(11c) 
(lid) 



Here, a <C A e ff is assumed. The partition function of the double-layer sine-Gordon model is thus identically rewritten 
in the form of a partition function for a gas of topological excitations, which we would like to call the "layered 
sine-Gordon gas" and which is given by 



v=o v J ' m=i J 1 J 



U2u — 1 ' 



( 2v 1 

cxp -b 2 £ 2 <Ja a f i Sr ' 



cri,...,tr 1/ = ±1 
<T„ +7 = -er 7 , 7 G{1,...I/} 



f A(r Q7 ) + (1 <5n„n 7 

)5(r Q7 )} 



(12) 



where the contact terms a — 7 are treated separately (the latter modification leads to a physically irrelevant renormal- 
ization of the partition function). The frequency b is inversely proportional to the temperature b 2 = 2ir/ (JcbT), and the 

Fourier amplitude y is related to the fugacity z by the relation z 2v /{v\) 2 = {y /2) 2v / {2v)\, i.e. y = 2z I ^ ^ J " J 
where (a) n = T(a + n)/T(a) is the Pochhammer symbol. 



C. Comparison of the layered sine— Gordon and layered vortex gas models 



The partition functions of the layered vortex gas and (Ti"2"|) of the layered sine-Gordon gas have the same 
structure. Therefore, the intra- and intcrlayer interaction potentials can thus be compared directly. A comparison of 
Eqs. (|5aj) and (jlOap reveals that for small distances (r <C A G ff), the intralayer potentials have the same logarithmic 
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behaviour for both models. This is not unexpected since in this case the vortices of a given layer are independent of 
the effects in the other layer. However, for large distances (r S> A e g), the intralayer potential is logarithmic for the 
gas of topological excitations of the layered sine-Gordon model in contrast to the the layered vortex gas, whereas the 
long-range intralayer potential is dominated by a linear term. The difference of the two models becomes even more 
significant if one compares the interlayer potentials which are different for the two models both in the short-range as 
well as the long-range regime (see Eqs. (|3b| and (|10bj) ). 

The significant differences of the long-range behaviour of the interlayer potentials strongly indicate different long- 
distance (infrared, IR) physics. We should note that the long-range behaviour of the potentials in Eqs. (|llb[) and 
(|1 ld[) generalizes to a leading asymptotics of the form — l/(2irN) ln(r Q7 /A e ff) for an TV-layer system with the interlayer 
interaction given in Eq. (|f 3[) below. Thus, the addition of more layers thus does not change the qualitative behaviour 
of the long-range potentials (linear versus logarithmic). We conclude that there is a strong indication that these 
models belong to different universality classes, and that the layered sine-Gordon model is not suitable to describe the 
vortex properties of Josephson coupled layered superconductors if a linear long-range potential between the topological 
defects is assumed. 



III. PHASE STRUCTURE OF THE TV-LAYER SINE-GORDON MODEL 



From a conceptual point of view, it is interesting to study the critical temperature as a function of the number of 
coupled layers, in the framework of an appropriate generalization of the double-layer model defined in Eq. ([5]) to TV 
layers. The discretization of the derivative term for the z-direction in the three-dimensional sine-Gordon Lagrangian 
results in a model of coupled 2D systems J33[, which has been called the TV-layer sine-Gordon model. It consists of 
TV coupled 2D sine-Gordon models of identical "frequency" b [13, [H| , each of which corresponds to a single layer 
described by the scalar fields (i = 1, 2, . . . , TV). Its bare action reads (see Eq. (2) of Ref. [22} ) 



s- 



NLSG 



f \l N N ~ 1 J N 

/ d2r 2 X^^) 2 + 2^ i+i ~ v?i - )2 +X^ cos ( 6< ^) 

i—1 i—1 i—1 



(13) 



We have implicitly defined the mass matrix m 



:JV 



of the TV-layer model, X^i* f (Vi+i 



action is invariant under a joint shift of all fields Lpi 



<Pi 



2ir/b applied to all layers 



= 1,2, 



IT 2 

2<p m N <£- 



The 



. , TV, a symmetry 



which corresponds to a single zero mass eigenvalue of the matrix rn 2 N - Indeed, after a suitable rotation of the mass 
matrix [12, , it becomes evident that the TV-layer sine-Gordon model consists of TV — 1 massive 2D and a single 
massless 2D sine-Gordon fields. The periodicity in the internal space spanned by the field is broken explicitly for 
the massive fields, and the spontaneous breaking of periodicity of the single massless mode accompanies the phase 
transition for small values of fugacities [13, H2] • 

The rotated TV-layer sine-Gordon model has already been investigated by the Wcgner-Houghton renormalization 
group method on the basis of the mass-corrected linearized scaling laws [33j and by a general perturbativc treat- 
ment [22j • Both approaches predict a linear increase of the critical parameter b\ with increasing number TV of the layers, 
according to the formula, b\ = 8ttTV, (see Eq. (35) of Ref. HU). Equation |(T2J) clearly implies b\ = 8ttTV = 27r/(fc B T), 
and we therefore obtain T c oc 1/TV for the TV- layer model. This decrease of the transition temperature is perfectly 
consistent with the general properties of the model in the limit of an infinite number of layers. Namely, one can 
intuitively assume that the single remaining zero-mass eigenvalue cannot make a decisive contribution to the phase 
structure of the model in the limit TV — > oo, with TV— 1 modes being massive. Indeed, in the limit of an infinite number 
of layers, one recovers the 3D sine-Gordon model which does not undergo any phase transition at all [30l. l32l. |33| . 



IV. SUMMARY AND CONCLUSIONS 



The main result of this paper (see Sec. [TTJ) is the indirect comparison of the layered sine-Gordon model to the layered 
Ginzburg-Landau theory: as we have shown, both models can be mapped to different gases of topological excitations. 
These are the layered vortex gas for the layered Ginzburg-Landau theory (see Eq. and the equivalent gas of 
topological excitations for the layered sine-Gordon model ("layered sine-Gordon gas, see Eq. (fT2"|) ). In general, we 
find that if a long-range confining linear potential is required for a description of the Josephson-coupled layered high-T c 
superconductors, then the system of coupled 2D sine-Gordon models is not suitable to describe the vortex properties 
of these materials: scalar-field propagators cannot provide linear potentials in two dimensions. For short distances, 
a logarithmic behaviour can of course be approximated quite well by a linear potential (see Ref. [13j |. ln(l + r) ~ r 
for r <C 1), but this observation is irrelevant for the phase-structure of a system, which is determined only by the 
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long-range interactions. In any case, we are led to the conjecture that the layered sine-Gordon and the layered vortex 
gas models belong to different universality classes. Using a renormalization group analysis of the generalized iV-layer 
sine-Gordon model as described in Sec. IIII| we find that the critical temperature of the layered sine-Gordon gas 
reads fulfills fee T c = (AN) -1 . This is inconsistent with high transition temperatures for multi-layer system and in 
strong disagreement with experiment [3^, HO, Hl[ . E.g., in Ref. [40| . for YBa2Cu307_,5 the single-layer (2D) transition 
temperature was determined as 30.1 K, and with N = 2 layers, the experimental result was Tktb = 58.2 K, suggesting 
T c (x N for a small number of layers. 

Let us conclude this paper with a perhaps somewhat surprising outlook. The decrease of the transition temperature 
with the number of layers is tied to the gradual "disappearance" of the "influence" of the only remaining zero-mass 
mode in the matrix of the Josephson-coupled layered sine-Gordon model, in comparison to the N — 1 massive 

modes, as N — > oo. If we choose the mass matrix differently, e.g., tp T M_ 2 ip = G (Y^ =1 a n (p„^j , with the (only) 

condition a;; = 1, then there are N — 1 massless modes and only one massive mode. In that case, we find (see [421]). 
T c oc ^fe^i and this result is in agreement with the analysis presented in Ref. [43| for magnetically coupled layered 
superconductors. In this case, the interaction potentials corresponding to Eq. ([TTj) between the topological defects 
have the same asymptotic behaviour as those given in Refs. for the magnetically coupled case. After all, a 

layered sine-Gordon type field theory with a suitable intcrlaycr interaction might prove to be useful for the description 
vortex dynamics in (magnetically coupled) layered systems, but not in the expected direction, which would have been 
the Josephson-coupled case. 
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